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Abstract 

In this article we present an analytic solution of the famous problem of diffraction and interference 
of electrons through one and two slits (for simplicity, only the one-dimensional case is considered) . 
In addition to exact formulas, we exhibit various approximations of the electron distribution which 
facilitate the interpretation of the results. Our derivation is based on the Feynman path integral 
formula and this work could therefore also serve as an interesting pedagogical introduction to 
Feynman's formulation of quantum mechanics for university students dealing with the foundations 
of quantum mechanics. 
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I. INTRODUCTION 



The quantum mechanical problem of diffraction and interference of massive particles is 
discussed, though without detailed formulas, by Feynman in his famous lecture notes.- A 
more exact treatment, though still lacking in detail, is in his book with Hibbs.- It was 
first observed experimentally by Jonsson in 1961. 3 Moreover, there are also experiments 
for neutrons diffraction for single and double slit (see^ and the references therein) and in 
quantum optics about interference between photons, (e.g.,-). The crucial point of this paper 
is to deal with different optical regimes: the usual Fraunhofer regime, and (less commonly 
taught to students) the Fresnel regime and intermediate regimes. Recall that these regimes 
depend on the distance between the slits and the screen, where the Fraunhofer regime 
corresponds to the case when the distance between the slits and the screen is infinite and 
the other regimes appear when this distance is finite, the intermediate and Fresnel regimes 
being distinguished by the value of the Fresnel number Np = 2a 2 /XL, where 2a is the width 
of the slit, L is the distance between the screen and the slit and A is the wave lenght of 
the electron. Thus, the purpose of this article is firstly to present the theory of the slit 
experiment using the Feynman path integral formulation of quantum mechanics, which may 
be of pedagogical interest as compared with the optical Young experiment (Section II, III, 
IV), secondly to give an analytical derivation of the final formulas for the intensity of the 
electron on the screen, and to analyse these formulas using some approximations based on 
the asymptotic behavior of the Fresnel functions^ occurring in these expressions (Section V). 
In particular we show how the physical parameters, especially the Fresnel number, affect the 
form of the diffraction and interference images. There exists some pedagogical papers about 
the multiple slit experiments (see e.g.- for experimental and numerical discussions and, 8 for 
theoretical discussions using path integral formulation), but surprisingly the approximations 
obtained in the Section V has never been published. 



II. FEYNMAN FORMULATION OF QUANTUM MECHANICS AND WHY IT 
COULD BE OF INTEREST FOR STUDENTS 

Students are often surprised to learn that under certain physical conditions the exper- 
imental behavior of matter is wave- like. In fact, this can present a didactical obstacle in 
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teaching the first course of quantum mechanics where the principle of wave-particle dual- 
ity can appear mysterious, especially with electron diffraction experiments for one- and two 
slits. Despite the interesting historical ramifications, students often have many metaphysical 
questions which are not answered satisfactorily in introductory quantum mechanics courses. 
A formal and complete solution of the electron diffraction problem based on Feynman's 
path integral approach may be helpful in this regard. It could demystify the diffraction 
experiment and clarify the principle of duality between wave and particle by analogy with 
wave optics. In addition, it could to be an interesting introduction to quantum mechanics 
via the Feynman approach based on the notion of path integral rather than the Schrodinger 
approach, highlighting the analogy between quantum mechanics and optics. Moreover, there 
is an interesting article that has recently been published in the European Journal of Physics 
Education that discusses using the Feynman path integral approach in secondary schools to 
teach the double slit experiment .— 

We begin by outlining Feynman's formulation of quantum mechanics. We recall some 
of the fundamental equations before studying the electron diffraction problem. For more 
detail and for historical remarks about this theory, see Ref.-,— ,— . We will give an equivalent 
formulation to that based on Schrodinger's equation, but which is related to the Lagrangian 
formulation of classical mechanics rather than the Hamiltonian formulation. Consider a 
particle of mass m under the influence of an external potential V [a; (£),£] where x(t) = 
(x^(t), x( 2 >(t), .., xW(t)) is the (d-dimensional, d > 1) coordinate of the particle at time t. 
The Lagrangian of the particle has the simple form 



where x(t) = dx(t)/dt denotes the velocity of the particle at time t. If the particle is at 
position Xi at time t« and at Xf at time tf > U then the action is given by 



In classical mechanics, the total variation of the action 5S for small variations of paths at 
each point of a trajectory Sx(t) is zero, which leads to the Euler-Lagrange equations well 
known to students of physics: 



L[x(t),x(t),t] = -x{tf -V[x(t),t] 




d dL[x(t),x(t),t) dL[x(t),x(t),t] 
Jt dx( a \t) dx&Ht) 



0, a — 1, .., d . 
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However, in quantum mechanics the Least Action Principle as described above is gener- 
ally not true. Thus, the concept of classical trajectory is also no longer valid: the position as 
a function of time is no longer determined in a precise way. Instead, in quantum mechanics 
the dynamics only determines the probability for a particle to arrive at a position Xf at 
time tf knowing it was at a given position Xi at time ti < tf. In other words, knowing the 
state of the particle at time ti, denoted by \xi, ti), the question is: what is the probability of 
transition from the state \xi,ti) to the final state \xf,tf)1 This probability is given by the 
square of the modulus of the so-called amplitude denoted (xf,tf\xi,U), i.e. | (x/, tf\xi, ti) | 2 . 
The expression for the amplitude in Feynman's formulation, though equivalent, differs from 
Schrodinger's firstly because the Lagrangian formulation is more general (there is not nec- 
essarily a Hamiltonian) and secondly because it does not refer to waves. In Feynman's 
approach the amplitude of transition is given by an 'integral' over all possible trajectories 
of a phase whose argument is the action divided by Planck's constant %, symbolically, 



where K(xf,tf;Xi,U) = (xf,tf\xi,ti) is the amplitude, where the action functional is 



where e = (tf — U)/n, x = Xi et x n = Xf. 

To understand this formula, we recall some ideas from Feynman's thesis^. Consider the 
particle at a point x at time t. Suppose that the particle changes position by an amount Sx 
during an infinitesimal time interval St. The action for this time interval can be written as 
S[x + 5, t + At; x, t] ~ L[x, %, t]8t. We define the amplitude of transition between the states 
\x, t) and \x + 5x, t + 5t) as 





and where j Dx(t) represent the path measure to be interpreted as follows: the integral 
formula is a short-hand for a limit of multiple integrals: 

^.*..*.^=n 1 f ^1 f ^=1 




(1) 




Xj Xj+1 



I i 



(Xo.to) tl 



(Xn.tn) 



ti ti+l 



H — I — I — I — I — h 

tn-l tn 



FIG. 1. Representation of paths between and (xf,tf), for the continuous case at the left 

(integral over continuum paths, see dU)) and for the discrete case at the right (integral over discrete 
paths before the continuum limit, see (JH)). in the left figure, the continuous line is the classical 
path, the dashed lines are two paths obtained by variations around the classical paths and the 
dotted line is an arbitrary path. In the picture at the right the discrete path for fixed n is shown 
together with the classical discrete path. 

Then, if we divide the time interval [£«,£/] into a sequence of small intervals = t < ti < 
t 2 < ... < £ n -i < tf = t n , where tk = ti + k(tf — ti)/n, the transition amplitudes are 
multiplied (because successive events are independent): 

(xf,tf\xi,U) = {xf,tf\x n -i,t n -i) ■ ■ ■ (x 2 ,t 2 \x 1 ,t 1 )(x 1 ,ti\xi,t i ) . 

Integrating over all possible values of the intermediate positions x^ leads to the formula (0Q), 
see Fig. [TJ The intuition associated with this formula is as follows: we integrate the phase 
over all possible trajectories of the particle; unlike the classical case, the particle can follow 
paths that differ from the classical trajectory which minimizes the action. (Paths far away 
from the classical path usually make negligible contributions.) This image, though based on 
the classical image of a trajectory, illustrates the change in the mathematical description of 
the particle (wave-like behaviour) which can no longer be represented as a material point as 
its trajectory is not clearly denned with certainty. 
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The integral equation describing the evolution of the wave function *&(x, t) (i.e. the state 
\x,t)) over time, is given by 

y(x',t')= / K(x',t';x,t)V(x,t)dx . (2) 

jR d 

It can be shown that the 'wave function' *&(x, t) is also the solution of Schrodinger's equation: 
— j-A x ty(x, t) + V(x, t)^f(x, t) = ihdt^f(x,t), where A x is the ^-dimensional Laplacian and 
d t = d/dt. For a proof of the equivalence between both formulations, see^,— . 

A useful example is the free particle amplitude calculation in one dimension. We simply 
replace the Lagrangian by the free Lagrangian, i.e. without potential, and use equation (pQ) 
to get : 

K (x', t'; x, t) = 1 e * . (3) 

17 y/2inh(t' - t)/m 

Let us derive formula (j3J) starting from the general formula (pQ) applied to a free particle in 

1-dimension. We divide the integral, for fixed n, we have to integrate the phase e iS, N,--,z„]/ft 

over the intermediate positions Xj, j = 1, .., n — 1 at times tj = je, e—(tf — tj)/n between 

x = Xi and x n = Xf, where S[x 0l ..,x n ] = multiplied by the constant 



rij=i ^/m/2iTrhej . We need to calculate, for a fixed number of subdivision n, the multiple 
integral: 

Kq^ (x n , t n ] Xq, to) 

I f°° dx\ f°° dx n ^i I . <^.m(xj — Xj-i) 



^2mTie/m J-oo y2infie/m J-oo \f2mTie/m ' y ~^ 2ne 



D m( \' H)i ) • (4) 



3=1 



We perform the integration one stage at a time (from j — 1 to j — n — 1), using the 
well known formula for the convolution of two Gaussians. Given two Gaussians g Gv [x) = 



-2 

) 



J- — e , v = 1,2, we have: 



g n (x - x')g a2 {x') = g y fis^s(x) . (5) 



Hence, in (jlj), we can calculate the integral with respect to x\, changing variable to x\ 
X\ — x , obtaining: 

1 f°° dX\ . m(^ 2 -^i) 2 ; m ( x i-xg) 2 

-e 2fie e 2he 



^2mhe/m J-oo ^/2mhe/m 



^2mh{2e)/ m 



_g 2fi(2e) 



Proceeding by recursion, we get (J4]) : 

^0 \ x ni t-nj -^0) L 0) — i — ir—. — e 1 ; > 

y 2mri(ne) / m 

and hence ([3]) in the limit n — > oo. 

Clearly the 3- dimensional propagator is a product of 1-dimensional propagators: 

K^ D \^', t'- it, t) = K Q {x', t'; x, t)K {y', t'; y, t)K (z', t'; z, t) 



( m \ 3 / 2 m \-£'--£ 
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\2mh{t'-t)) ' lJ 

where x = (x, is the position vector in 3-dimensions. Note that for reasons to be 
discussed below, we shall not make use of the 3-dimensional propagator in the example of 
election diffraction for one and two slits. We remark for completeness, that the propagator 
(J3]) can be calculate by other methods, for example.— 



III. APPLICATION TO THE PROBLEM OF ELECTRON DIFFRACTION AND 
INTERFERENCE 

Consider an electron source at (x, y, z) = (0, 0, 0), and two slits at z = D of width 2a and 
centered respectively at x = +b and x = —b, c.f. Fig. |2]. For the diffraction experiment with 
a single slit, just replace the system with a slit of size 2a centered at x — 0. At z = D + L is 
a screen on which electrons are recorded (or some other recording device). For more details 
about the experimental realization of the system, see Ref.- Note that we neglect gravity for 
this problem. In addition, we assume that in the direction orthogonal to the plane in Fig. [2] 
(the x — z-plane), the slot is long enough to neglect diffraction effects; we consider only the 
horizontal (in-plane) deflection of the beam in order not to complicate the formulas. 

Then we can reduce the dimension of the propagator ([6]) by integrating over y: 

/oo 
dyK^ D \^\t'^,t) 
-oo 

m - ml^'-^l 2 

e 2E(f-t) ) (7) 



2i7rh(t' - 1) 

where "r^ = (x, z) is the position vector in the two dimensional plane orthogonal to the 
y-axis. 

One can suggest different models for the slits given by distribution functions (e.g. Gaus- 
sian functions or 'door' functions). We focus on the more realistic model of 'door' functions 



(for the Gaussian function model, see Ref.-): 

{0 w > b + a, w < b — a 
( 8 ) 
1 b — a < w < b + a 

The question is: what is the probability of finding the electron at the point x on the screen 
knowing that it started at the point (x — 0, z — 0)? More precisely, suppose that the source 
emits electrons in large numbers, though small enough so that the distance between elec- 
trons is such that interactions can be neglected (no correlations). What is the intensity of 
electrons on the screen as a function of position xl The two questions are related since under 
these circumstances the motion of the electrons is independent (without mutual interaction, 
given that the density of the beam is very low). Indeed, the intensity curve is obtained by 
simply multiplying the probability curve for an electron by the number of electrons emitted 
per unit time. As explained above formula ([TJ, the probability is given by the square of the 
amplitude which we will compute using the propagator flSJ). 

As explained in (J7J), we now take d = 2. In fact, we now argue that we can reduce 
the dimension further. For the ^-direction, we really onght to use the two-dimensional 
propagator but we can reasonnably consider (as it is usually done^,-) that the propagator is 
a product of two independant propagator in both directions x and z, where the last one is 
equal to a constant. 

To see that, let us discuss the conceptual experiment for one slit.- We consider that 
the problem is divided in two separate motions, one starting from the source to the slit 
during the time T and the other one starting from the slit to the screen during the time r. 
Then, we would like to compute the probability amplitude for the electron from the source, 
at the intial position (x — 0, z — 0) at the time t = 0, to go to the screen, at the final 
position (x,z = D + L) at the time T + r, knowing that it goes through the slit, at the 
intermediate position (w, z = D), b — a < w < b + a, at the time t = T. In fact, by the 
laws of quantum mechanics, there is no reasons to separate the motions in two independant 
parts, since we don't know where is the particle at the time T.In other words we don't know 
when the particle goes through the slit. Nevertheless, we can consider that this classical 
image is appropriated to study the problem. Indeed, the electron have in the ^-direction 
the momentum p z = hk z (k z is the wave vector), which is related to the classical velocity 
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v z = D/T = L/t, where D is supposed to be very large compared to the dimensions in 
the x-direction, more precisely, x,a,b <C D, L. In addition we suppose that the wave lenght 
A, which is approximatively equal to the ^-direction wave lenght A ~ X z = 2irh/(mv z ), is 
small compared to the distances A <C D,L. Thus, the motion is approximatively classical 
in the z-direction and we can separate the problem as two independant motions. Notice 
that, quantum-mechanically, it is possible for the particle to go through the slit severals 
time before strike the screen,— but that the probability of this event is relatively small. 

Now, let us compute the amplitude of the transition for the particle starting at the point 
(x = 0, z = 0) at the time t = 0, going through one slit at the position (w,z = D), b — a < 
w < b + a at the time t = T and arriving at the position (x, z = L + D) at the time t = T + r: 

/oo 
-oo 

K^ D) ((X, L + D),T + r; (w, D), r) K^ 2D) ((w, D), r; (0, 0), 0) . 
Hence, the explicite formula is given by: 
K { ™\{x,L + D),T + t; (0,0), 0) = 

b+a 



e 2RT e 2Bt / g Ihr e 2BT 

, = , == / dw = = . (9) 

\J1mKr jm ^/ImYiT jm Jb-a y 'ImTvr / 'm ylmKT jm 

Consequently, the two-dimensional propagator is the product of two independant one- 
dimensional propagators in the x and in the ^-directions, and since the propagator in z- 
direction is a constante (see the right hand sides of Qj): 

K (Zy D T '0 0) 2RT ( ) 

sJlrnKT jm ^/HtxTit jm 

we can reduce in the one- dimension's x-direction the problem: 

rb+a ■ m(x-w) 2 ; mw 2 

K x (x,T + r;0,0)= dw r^-^ ^=^f= • (") 

Jb-a y liixar jm ^ztTthl jm 

Then we can reduce the dimension of the problem keeping only the x-axis propagator, as 
you can see in the References^,-. In the Reference,- they consider the slits in the x-y-plane, 
in two-dimensions and not in one-dimension as we do in the present article. 
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If we take the limit a — > oo (infinite slit), by (JSJ) and by ([5]) we get: 

j mD 2 - ml? ,- mp-n)) 2 -mii£ 

^(2D)^ ^ y '(0 0) 0) e 2ST e 2Ht d 2hT e 26T 

y^2mhT/m ^2mUT jm J-oo \/2iixfiT jm ^2mTiT / 



■ mL 



^j2inhT /m ^j2iix%T /m ^j2iixh(T + r/m) 
Then we find the amplitude corresponding to the transition between the position (x — 0, z — 
0) at the time t = and the position (x, z = L + D) at the time t = T + r, knowing that 
the position in the z-direction is D at the time t = T : 

Kg D \(x,L + D),T + T;(0,0),0) 

= K 0>x (x, T + r; 0, 0) K 0> * (L + D, T + r; £>, r) K , 2 (A T; 0, 0) , 

where K 0tX (X', X, t) = (m/2iixh(t' - i)) 1 / 2 exp (im(X' - X) 2 /(2h(t' - t))), X = x,z, are 
the 1-dimensional free propagators. Then, integrating over the intermediate positions D 
between — oo and +oo, we find the two dimensional free propagator ([7]) between (x = 0, z — 
0) at t — and (x, z = L + D) at the time t = r + T: 

/oo 
dD K (2D \(x,z = L + D),T + t; (0,0), 0) 
-oo 



oo 



■ mD i ■ m(z-D) 



2 



e % ^2nr~ e — 2HT — e 2fi ( T +^) 

oo \/2mhT lm \/2inhT /m \/2iirh(T + r)/m 
rn 



^j2m%T /m ^j2in%T jm ^j2iitfi(T + t)/i 

■ m(x 2 +z 2 ) 
e l 2H(T+t) 



= K^ D \(x,z),T + t; (0,0), 0) , 
to get the result, we have used the convolution formula between two Gaussians §5§. 



IV. THE EXACT RESULT IN TERMS OF A FRESNEL INTEGRAL 

Now we want to compute the amplitudes Ai(x; a,b), A2(x;a,b) at each point x on the 
screen, using the Feynman formulation, and then to add both amplitudes to get the total 
amplitude A(x; a, b) and finally to take the square of the modulus, obtaining the probability 
P(x;a,b) = \A(x;a,b)\ 2 . 

By ffTT|) . the formal expression for A\(x; a, b) is : 

/+oo 
dw X[b-a,b+a](w)K {x,T + r;w,T) K (w,T; 0,0) , (12) 
-oo 
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FIG. 2. Schematic depiction of the apparatus with the source, the two-slits and the screen. We 
have represented four differents paths going through the slits, from the source to the screen. The 
dashed lines represent three classical paths and the thin dashed line a curved path. The Feynman 
path Integral formulation allows one to compute the propagator summing over all different paths, 
see equations CD), ©, flE}, (US]), jEJ), (ITS]) . 

where T is the travel time of the electron from the source to the slit and r from the slit to the 
screen. Recall that to obtain formula ( TT2|) . we used a similar argument to that which enabled 
us to write the formula (JT]), writing the integral over all possible paths as the product of 
independent amplitudes at successive times. However, in this case, at time T, we have to 
integrate over a finite interval (the slit), which results in the more complicated expression 
( I12p rather than a Gaussian: 



Ai(a;;a, &)= / rfw e 2 B r =e aar . (13) 

Jb-a y2m%T /m ^ImKT jm 



Notice that : 



mix — w) 2 mw 2 ,171 m w a; N0 ma; 2 

-(t^z + ^f (y- 1 i ,J + 



2frr 2/iT v 2ftr 2frT yvy l + r/T y 2ft (T + r) 
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and hence, 

■ mx^ Li/ 

e'a(T+T) f b + a t + t ( T + T X 

Ai(x;a, b) = — ; / dw* , — -. — exp I i{— — — -. — )(w — 



2mh(T + r) /m J b - a \ ZmhTr/m * \ 2hTT/m 1 + r/T 

mx^ (1) / \ 

e l 2h(T+r) /■«+ W 



^(2i)Hh(T + r)/m J a w {x) 

where 



dw'ex P (—w' 2 ), (14) 



a±(x;a,b) = \l ^ + T j- (b ± a) = J — ^— . (15) 

In ( fT4l) . we see that we have the integral of a Gaussian with complex argument. Decomposing 
the integral in real and imaginary parts, we get two integrals of cosine and sine functions 
respectively, with second degree polynomial arguments. These integrals are the well-known 
Fresnel functions^: 



C[u] = dw cos 
Jo 

S[u] = dw sin ( 
Jo 



,nw 2 « 



TTW 



2 

2 



'0 

Thus we obtain explicit analytical expressions for the amplitudes: 



g'2fi(r+T) 

Ai(x:a,b)= — = x 

v /(2i)Hh(T + r)/m 

(C[a + (x; a, b)] — C[a_(x; a, b)] + iS[a + (x; a, b)] — iS[a-(x; a, b)]) , (16) 

A 2 (x; a, b) = A\(x\ a, — b) . (17) 
For two slits, we can compute the total amplitude by summing the amplitudes for both slits: 

A(x; a, b) = Ai(x; a, b) + A 2 (x; a, b) . (18) 



V. PHYSICAL PARAMETERS, APPROXIMATIONS AND INTERPRETATIONS 

Recall that the length of the slits and the distance between them is small compared to 
the horizontal distances D and L, hence we can assume that the wave length of the electron 
A = h/mv is approximatively given by h/mv z , where v z = L/t = D/T. We will use this 
expression for the wave length in the following. 
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A. Diffraction by a single slit 



We can write for the single-slit case (of size 2a), the analogue of the function defined by 

a(x; a) = ^/N F (a)^l + L/D (l - f 1 + (19) 

where Np(a) = 2a 2 /XL is the Fresnel number. 

We can now easily compute the single-slit diffraction probability: 

p( 1SHt \x;a) = \A 1 (x;a,b = 0)\ 2 
= 2X ^ + D) ([C{a{x- a)) + C(a(x; -a))] 2 + [S(a(x; a)) + S(a(x; -a))] 2 ) . (20) 

Let us introduce the following parameters rj = 1 + L/D and 7 = 77 — 1. We can then plot 
the functions (120]) for various values of these parameters, see Fig 3. Note that the Fresnel 
functions in ( J201) behave differently depending on the value of the Fresnel number Np(a), 
esp. depending on whether it is greater or less than unity. To understand these differences 
explicitly, we will analyse the asymptotic behavior of these functions for different regimes of 
Np(a) using the known asymptotics of the Fresnel functions: (see^) 

C(±u) ~ ±- + — sm — , u > 1 , 

2 7TU 2 

S(±u)~ ±--— cos— , w>l. (21) 

2 7TU 2 



If we have: 



x 1 
- 1 > 



ar) ' ^N F {a)r) ' 
we get the asymptotic behavior for the functions defined by (TH)]) : 

a(x; a) <C —1 and a(x; —a) ^> 1 ±ct(x; =Fa) 3> 1 , (22) 

and so by (I2T]) and (|22|) . the asymptotics formulas of the Fresnel functions in (1201 are given 
by: 

_ r , , 1 1 . ,7ra;(2;; ±a) 2 , 

Cb i; ±a) — ±- H ; r sm ( K — '-), 

L V ' n 2 na(x;±a) K 2 h 

„x , , 1 1 . ,ita(x; ±a) 2 . . 

S[a (* ;±a )]. ± --_^ sm (^^). (23) 
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Then we get: 

_, r . m ^.r / m 1 • ^oi(x; a) 2 . 1 . .-Ka(x;-a) 2 . 

C ax; +a + 6 a x: —a ~ ; sin H ; sm 

L v ' yj L v ' yj vra(x;a) v 2 ' vra(x; -a) v 2 



or / m nt i m _1 ,7ra(x;a) 2 . 1 ,7ra(x; -a) 2 

5 a x; +a + 5 a x; -a ~ — r cos 1 , cos {—^. '-) 24 

7ra(x;a) 2 7ra(x; — a) 2 



Applying the Fresnel function asymptotic forms ( |24"1) to ( 1201) and using the definition ( |T9f ). 
we deduce than if N F (a) <C 1 and if (x — arfj/ar} 3> 1/ y/Npiajr) ^> x — ar? ^> a/ XL/2, we 
get the following asymptotic formula: 

p( 15M )(x;a) ~ 4^ f -r^ + ^^sin 2 (7riV F (a)-) 1 , (25) 

Moreover, since Np(a) <C 1 and so ^- 3> 1, we have another asymptotic form (large 



distance on the screen): 



p( 1SHt \x; a)^^L sin 2 (vriV F (a)-), (JV F (a) < 1, x/ar, > 1) , (26) 

In this case we are in the so-called Fraunhofer regime analogous to plane wave diffraction in 
optics.— In fact, notice that the distance between fringes is a/Np(a) = \ z L/2a, c.f. Fig 3a. 

Both approximations ([25]) and (126]) are also valid if Np(a) is of the order of unity (this is 
the intermediate regime) provided that x ^> arj. This means that the pattern on the screen 
far from the position of the first lobe is well approximated by equations (1231) and (1261) . see 
Fig. 3b. 

On the contrary, if Np(a) ^> 1, we get different asymptotics given by: 
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P (ism )(x) ^ 7 vM sin(|iV F (aMl-^) 2 ) sin (fA^Ml + 
1 ' V\ 2a 2 7 r v ^(a - 2) 2vr v ^(a + 2) 

2 / Vj\gg cos(fiV F (aMl-^) 2 ) cos (f ^(0)77(1 + ^|) 2 ) \ 2 
%V 2a 27r^(a-|) 27^ + 1) J ' < ^ W 

^(*)- ^ + ^ - 2 (^(«#) , N > ^ . (28) 

since the asynptotic behavior of the functions ( 1T9|) are: 

a(x; ±a) > 1, if iVp(a) > 1 and |x| < ar/ , (29) 
±a(x; =Fa) ^> 1, if Nf(o) ^> 1 and |x| > ai] , (30) 
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so if \x\ < ar], by (129]) and (l2Tj) . we get: 

_, r . , 1 1 . ,-naix] ±a) 2 
C[a(s; ±a)] c - + -^—^ sm ( ), 

then by ( 129]) and (l3Tj) . we have: 

„ r , „ r , X1 1 . ,7ra;(x;a) 2 . 1 . .ira(x; —a) 2 . 

C[a(x; +a + C[a(x; -a ~ 1 + — sm 1 + — sm (— '-), 

7ia{x;a) 2 7ra{x; — a) 2 

or / m or , m 1 .7r«(x;a) 2 . 1 .7ra(x;-a) 2 . 

i> +a) + 6 afx; -a) ~ 1 7 cos ( ) -. cos ( ) , 

ira(x;a) 2 7ia(x; —a) 2 

then we get ( 1271) . 

If |x| > ar], by ( 130|) and ( I2TI) . we get the same approximations as (|23l) and ff2^|) . then we 
obtain fl28|) . 

Note that the function ( 1271) oscillates rapidly in the interval [—a, +a] (esp. near the edges) 
around a constant value Np(a)'~f /2a 2 r] = 1/(X Z (L + D)), whereas for \x\ > ar], the function 
( 128|) decreases rapidly to 0. Hence p( 1Sht ) tends at large N F (a) to the 'door' function defined 
by (JHJ), as might have been expected, see Fig. 3c. 

B. Comment about the probability interpretation 

We can see that (120|) has the physical dimension of the inverse of a length squared 
and so it is neither a probability nor a probability density. This apparent problem can, 
however, be seen to be a matter of interpretation by looking at the formula ()2]). Indeed, the 
probability density for the diffraction problem is given by \^{x, T + r) | 2 , where ^(x, T + r) 
is a normalized wave function, i.e. J R1 dx \^(x,T + r)| 2 = 1. Therefore, we must choose an 
initial wave function (at time t = 0) which is also normalized so that its square modulus 
describes the probability distribution of the electron in the Xq plane. Essentially, what we 
have done in ( 120]) is to take the initial wave function to be a delta-function, whereas it should 
be the initial probability distribution which is a delta-function, i.e. ty(x, 0) should be 'the 
square root of a delta-function'. 

To make this clearer, consider for example a wave function at time t = given by 
the square root of a Gaussian <j) a (x ) = ( 27ro -2)i/4 e ~ i ° ? ? so t na t f Rl dxo \4>a{xo) \ 2 = 1 and 
|0o-(^o)| 2 — > $( x o) as o~ — y 0. In this way the wave function obtained in (134|) is properly 
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FIG. 3. Diffraction curve for a single-slit, with r\ = 2. The abscissae are the distances in units of a 
and the ordinates are the relative populations. We have Np(a) = 0.01 for the Figure (3a), 0.5 for 
(3b), and 100 for (3c). 



normalized so as to get the probability of the presence of the electron at the point x on the 
screen by taking the square of the modulus. Indeed, the wave function at time t = T, i.e. 
at the position of the slits, is given by: 

(f) a (x,T) = / dx K (x,T;x ,0)(j) a (x ) , (32) 
Jr 1 

where K (x,T + r;xo,0) is the free propagator defined by the equation ([3l). Using the 
identity 

/ K (x,T;x ,0)*K (x,T;x' ,0)dx = 5{x - x' Q ) (33) 
Jr 1 

we have that <f) a (x,T) remains normalized, i.e. J \4> a (x, T)\ 2 dx = 1. 
The wave function at time t = T + r is given by: 

V a (x,T + T)= [ dyK^x.T + r-y.T^y.T). (34) 

J — a 
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Now, the quantity of interest is the conditional probability (density) for the electron to be 
at the point x on the screen at the time T + r given that it was in the interval [—a, +a] at 
time T, i.e. given that it passed through the slit: 

l\P (x T + t)\ 2 

P a (x, T + r | y G [-a, +a),T) ee ' j ' } » . (35) 

J - a dy\My, T )\ 2 



after which we wish to take the limit a — > 0. Using the relation (1331) one can see that the 
condition probability ( |35l) is normalized so that this procedure just amounts to division by 
a normalization factor. Thus 

/ dx\y a {x,T + T)\ 2 = [ + dy\My,T)\ 2 . (36) 

JRI J-a 

'2 

To take the limit a — > 0, note that (p a (x ) = (8no- 2 ) 1 / 4: g ay ^(x ) where ^^(^o) = ^/j— r * s 
a normalized Gaussian of variance ay/2, and hence ^^(xo) — > ${ x o), °~ ~ > 0. Multiplying 
top and bottom of AMI) by (87TCT 2 ) 1 / 2 we thus have 

lim CT ^ P CT (x G A, T + r|y G [-a, +a],T) = 
(Stt^-^I^^^ + t)! 2 



lim 



>o (8^-1/2 /_^dy|^(y,T)|2 

1 2 

I / dx K(x, T + r; g , O^^Xp) | 

/-a d f |/rfxo-^o(y,T;x , 0)5-^^2(^0) I 
|if(x,T + T;0,0)| 2 
~ dy\K (y,T; 0,0)\ 2 

= ^P 1Sht (x;a) } (37) 
K(x,T + r;0,0) = / dyK (x, T + r; y, T)K Q (y, T; 0, 0) (38) 

./ — a 

is the propagator through the slit and P^ 1Sht \x;a) = \K(x,T + r;0,0)| 2 is given by (1201) . 
Note that this now has the correct dimension of an inverse length. 



where 



C. Interference and diffraction for two slits 



Similarly, one can find the two-slit diffraction probability formula using ( JT6l) . (JIT]) : 

P {2Slit) (x;a,b) = P 1 (x;a,b)+P 2 (x;a,b) + I 12 (x; a,b) , (39) 
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with the diffraction terms 



Pi(x;a,b) = \Ai(x; a, b)\ 2 
= [\C(a + (x] a, &)) — C(a_(x; a, 6))] 2 + [5(a+(a;; a, 6)) — 5(a_(a:; a, &))] 2 ) , 

P 2 (x;a,b) = \A 2 {x;a,b)\ 2 = Pi(x;a, -b) , (40) 

and the interference term : 

Ii 2 (x; a, b) = Ai(x] a, b)A 2 (x; a, b)* + A 2 (x; a, b)A\{x; a, b)* 

= (\C(a + (x;a,b)) — C(a^(x; a, b))][C(a + (x; a, — b)) — C(a^(x; a, — b))\ 

\Lr\ 

+ [S(a + (x; a, b)) — 5*(a_(x; a, b))][S(a + (x; a, —b)) — S(a^(x; a, —&))]) • (41) 

Notice that there is an additional term compared to the single-slit case called the interference 
term, which is of course quite similar to that in optics. 14 This results in a modulation effect 
of the curve given by ( 1391) by the sum of the diffraction terms ( [401) (modulo a multiplicative 
factor), see Fig. 4. 

Let us define the Fresnel numbers 

N F (a) = 2a 2 /X z L, N F (b) = 2b 2 /X z L and N F = 2ab/X z L = y/N F (a)N F (b)/2. 

Assume that the distance between the slits is large compared to the size of the slits b^> a. 
In the experiment considered one fixes both parameters a and b and varies the distance 
between the screen and the slits (keeping the same value for 77). Notice that because b 3> a, 
N F 3> 1 does not necessarily imply that N F (a) 3> 1. Thus we will see that both parameters 
play different roles. 

First, we establish the asymptotics of (J3HD for different asymptotic value of N F (a). Under 
the condition N F (a) <C 1 and at large scales \x — br]\ 3> at] et \x + brj\ ^> arj, we get 
similar expressions for Pi (a;; a, b) and P 2 (x; a, b) as (!26l) . We have to compute the asymptotic 
expression for the interference term. This yields 

p(23 ""^ 6 > - ^Bw sm2 ~ £>) + ^rw 8in2 i M + 

- 2/ 2 7 ,2 2^ (cos(27r(JV F + JV F (o))-) - cos (27^(1 + —))] 
7T z [x z — b z 7] z ) \ a arj J 

~ ^1^2) ( cos (M^-^(o))^)- cos (2^(1-^)))) . (42) 
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One can observe that there are two phases: the separated phase (Np 3> 1) and the mixed 
phase (Np < 1). At the same time there are the Fresnel and Fraunhofer regimes depending 
on the values of Np(a) as explained above. The distinction between two phases is purely 
geometric and characterizes the separation respectively mixture of diffraction curves. Indeed, 
similar to optics, to observe the two diffraction curves separately, the fringe modulation 
X z L/2a must be less than the distance between the origins of the two curves (being centered 
in ±br]) because otherwise both curves are mixed. Thus the criterion is written X z L/2a < brj 
and therefore Npi] > 1. This obviously unlike the case Np <C 1 where the two curves are 
combined (one added to the other) and where we observe modulation interference. 

To see this more formally, consider firstly the case Nprj <C 1. If \x\ > XL /2a then \x\ 3> br\ 
and we can give an approximation of (j42p . Indeed, the first two terms are approximately 
equal and contribute 

■ sin (——£). 



H 2 X 2 XL 

In the last two terms we develop the cosine functions and get 

2nb 2ixa 2irabr] 2nb 

cos (——x) cos (— -x) — cos ( — — — ) cos (——x) 
XL XL XL XL 

,2ixb 2 2ixa 
« — 4 cos(— — x) sin {——X), 
XL XL 

where we used cos (^rf 2 ) — 1 because Npt] <C 1. Adding the terms we obtain 



P( 2sht \x; a, b) ~ Vr sin cos T7 1 ■ 4 3 

n z x z XL XL 



This is the familiar optical formula: see Reference,— formula (10) Chap. VIII-6. It shows 
that the diffraction curves are modulated by interference fringes. The distance between two 
interference fringes is of the order of X z L/2b whereas that between minima of the diffraction 
curves is of the order of X z L/2a ^> X z L/2b, see Fig. 4a, Fig. 4b (far from the first lobe) and 
Fig. 4c. 

Secondly, consider the case that Npi] ^> 1 (while still assuming Np(a) < 1). If x — br) > 
—XL/2a (or x + bq < XL/2a) then (x — br\)~ 2 (x + brf)~ 2 (respectively (x — brf)^ 2 
(x + br])' 2 ) so one of the two terms is negligible in the respective domain. Moreover, in both 
cases, the interference term is small compared to the diffraction term since (x — br])' 2 ^> 
(x 2 — b 2 !] 2 )' 1 . The total probability is therefore approximatively equal to a sum of the 
two diffraction curves centered at ±677 modulated by an interference term which oscillates 
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rapidly with a relatively small amplitude, c.f. Fig 4c : 

P^\x;a,b) ~ 2 J sin 2 f vriV^l T f )) + Q( 2 1 2 ), |*±H > ^ (44) 
7r z (x^f bri) z \ or] J x A — b z rj z 2a 

Now, consider N F (a) ^> 1. Since, unlike the previous cases, we do not need special 
conditions for the position x on the screen, we find similar formulas to (I27p and (I28p for the 
direct terms Pi and P2, except that x is replaced by x — brj for the slit centered at +b and 
by x + brj for the slit centered at —b. 

For the interference term, inserting the asymptotics (I2ip into (14"T!) above results in the 
sum of two terms, one being the product of differences of cosin-functions, the other the 
product of differences of sin-functions. The problem is obviously symmetric about x = 0, so 
we need only consider the case x > 0. Then there are again two cases: 
(i) \x — b\ > a; and (ii) \x — b\ < a. 

In the first case, both terms decrease like l/{x + b)(x — b) with various fluctuating factors as 
in fT42|) . In the second case, Pi behaves as in ff271) but centred around x = b, and the other 
terms are negligible. We do not write the asymptotic formulas explicitly because the result 
is simply the observation that in this case we obtain a sum of two separated diffraction 
curves in the Fresnel regimes, i.e. curves that tend to the door functions in the limit, see 
Fig. 4d. 

Notice an interesting behavior of the interference pattern in Fig. 4c, where we see that the 
interference amplitudes are very small compared to the diffraction amplitude inside a band 
25^|x/a|^75, so that there are no interference fringes. This is also discernible in Fig. 4b 
of the Reference^ which corresponds to the calculated two-slit diffraction images, where one 
can observe the absence of fringes in a band. However, this phenomenon is not apparent on 
the corresponding experimental image. This is probably due to the difference in defocussing 
between the calculated and experimental images, see Fig. 3b and Fig. 4b of the Reference.- 
Indeed, the existence of such a band is quite sensitive to the value of the parameter Np(a). 

VI. CONCLUSION AND REMARKS 

- We have briefly presented the Feynman approach to quantum mechanics, based on the 
Lagrangian formulation of classical mechanics, and the associated change in paradigm in 
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(4a) (4b) 




FIG. 4. Interference curves (I39p from two slits, with b/a = 13, 77 = 2. The abscissae are the 
distances in units of a and the ordinates are the relative populations. We have Np(a) = 0.001 for 
the Figure (4a) , 0.015 for (4b), 0.12 for (4c) and 6 for (4d). 



the transition from classical to quantum mechanics. We note that in this approach, the 
transition from classical to quantum mechanics is quite natural because it relies mostly on 
concepts well known to students of analytical mechanics and does not confuse particle and 
wave behavior. It thus avoids some metaphysical questions and leads directly to the solution 
of the diffraction and interference problems above, and hence to a better understanding of 
the quantum mechanics of such quintessential phenomena. This justifies introducing the 
Feynman formulation at an early stage especially as the semi-classical approach that is 
often used in a first course relies on the idea of quantification of the action. A parallel 
introduction of the Feynman integral thus makes sense as it clarifies the passage classical to 
quantum. 

- Secondly, it seemed of interest to derive explicit formulas for the problem of diffraction 
/ interference by one or two slits, and to discuss the results based on the physical parameters 
of the system, notably the Fresnel numbers and the distance scale at which we observe on 
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the screen. The properties of the diffraction and interference patterns are not apparent from 
the exact formulas (1201) fl39l) . so it is useful to establish asymptotic forms f l25|) ff26|) fl27|) for 
the case of one slit, and ( l4"2"j) for the case of two slits. 

We summarize the various conclusions. In case of a single slit: 

- If Np(a) <C 1, this is the Fraunhofer regime for which the distribution curve is similar 
to the plane wave case, c.f. equations f )25|) . ff26|) and Fig. 3a. 

- If Np(a) ^> 1, this is the Fresnel regime for which the diffraction curve approximates 
the form of the slit, c.f. equations ( 1271) . ( 1281) and Fig. 3c. 

- If Np(a) ~ 1 one is in the intermediate regime for which there is a spreading around 
the center of the electronic distribution and we find the case of Fraunhofer distances on the 
screen, c.f. equation (120]) and Fig. 3b. 

In the case of two slits of width 2a, and separated by a distance 2b with b ^> a, we 
can make similar distinctions as in the one-slit diffraction case but there is also a transition 
between two phases dependent on the optical resolution: 

- If Np <C 1, one is in the mixed phase, i.e. we observe an interference curve modulated 
by a diffraction curve for a slit of size a in this case Np(a) <C 1, then we are in the regime 
of Fresnel, c.f. equations ( l4"2l) . ( 143]) and Fig. 4a. 

- If Np ^> 1, one is in the separated phase, and there are two interference curves (the 
interference amplitudes are lower) modulated by the diffraction curves corresponding to both 
slits, each curve being centered respectively at ±677; the shapes of the diffraction curves 
modulating the signals of each of the slits depends on Np(a) and are similar to the case of a 
single slit as summarized above (with three regimes: Fresnel, Fraunhofer and intermediate), 
c.f. (SID, Fi §- 4c for n f(o) < 1 and Fig. 4d for N F (a) > 1. 

- If Np ~ 1 we observe a separation between two interference curves, modulated by the 
diffraction curve corresponding to one slit at the intermediate regime, see Fig. 4b and Eq. 

Note that the fringes corresponding to the diffraction are at a distance XL /2a and those 
for interference at about XL /2b. The analytical properties of our asymptotics of two slits 
do not permit us to estimate these distances more exactly, but by analogy with optics they 
may be considered adequate. 
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In perspective, I suggest to take into account the quantum-mechanical way in the in- 
direction to solve the problem completely. Indeed, recall that as we discussed in the Section 
III, we consider in this article that the problem is separated in two motions, one between the 
source and the slits and the other one between the slits and the screen, which is rigorously 
not true. This is a challenging task since we have also to compute the loop path corrections,— 
but it could be an interesting contribution to the European Journal of Physics. 

ACKNOWLEDGMENTS 

I would like to thank the Professor Tony Dorlas for discussions, encouragements and 
English corrections of the manuscript. 

REFERENCES 



* |mbeau@stp.dias.ie| 

1 R. P. Feynman, R. B. Leighton, and M. L. Sands , The Feynman Lectures on Physics (Addison- 
Wesley, Reading, MA, 1963). 

2 R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals (New York: McGraw- 
Hill), 3rd. ed. (1965) 

3 C. Jonsson, "Elektroneninterferenzen an mehreren kiinstlich hergestellten Feinspalten", Z. Phys. 
161 (4), 454-474 (1961). 

4 Anton Zeilinger et al, "Single- and double-slit diffraction of neutrons", Rev. Mod. Phys. 60, 
1067-1073 (1988) 

5 R. L. Pfleegor and L. Mandel, "Interference of Independent Photon Beams", Phys. Rev 159 (5), 
1084-1088 (1967). 

6 M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, (Dover, New York, 
1965). 

7 Stefano Frabboni et al, "Two and three slit electron interference and diffraction experiments," 
Am. Journ. Phys. 79 (6), 615-618 (2011). 

23 



A. O. Barut and S. Basri, "Path integrals and quantum interference", Am. J. Phys. 60, 896-899 
(1992). 

9 M. Fanaro, M. Arlego, M. R. OTERO, "A Didactic Proposed for Teaching the Concepts of 
Electrons and Light in Secondary School Using Feynmans Path Sum Method", European J of 
Physics Education 3:2, 1-11 (2012) 

R. P. Feynman, Feynman's Thesis: A New Approach to Quantum Theory, edited by Laurie M. 
Brown (World Scientific, Singapore, 2005). 

1 R. P. Feynman, "Space-time approach to non-relativistic quantum mechanics," Rev. Mod. Phys. 
20, 367-387 (1948). 

2 H. Kleinert, Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial 
Markets (World Scientific Publishing Co., Singapore, 4th edition, pp. 1- 1547, 2006). 

3 H. Yabuki, "Feynman Path Integrals in the Young Double-Slit Experiment", Int. J. Theor. Phys 
25:2, 159-174 (1986) 

4 M. Born and E. Wolf, Principles of Optics: Electromagnetic Theory of Propagation, Interference 
and Diffraction of Light, 4th ed. (Pergamon, Oxford, 1969). 



24 



